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absolute age indicator from the presence/absence of red 
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ABSTRACT 

The determination of age is a critical component in the study of a population of stellar 
clusters. In this letter we present a new absolute age indicator for young massive star 
clusters based on J-H colour. This novel method identifies clusters as older or younger 
than 5.7 ± 0.8 Myr based on the appearance of the first population of red supcrgiant 
stars. We test the technique on the stellar cluster population of the nearby spiral 
galaxy, M83, finding good agreement with the theoretical predictions. The localisation 
of this technique to the near-IR promises that it may be used well into the future with 
space— and ground— based missions optimised for near-IR observations. 
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1 INTRODUCTION 

In order to study a population of stellar clusters, i.e. to de- 
rive cluster dispersal times, relate cluster to star-formation 
events in the galaxy, or estimate their impact on the sur- 
rounding interstellar medium, the ages of individual clusters 
must be accurately measured. For this purpose, there have 
been a variety of techniques proposed in the literature to 
estimate the age of young massive clusters (YMCs), both 
in the Galaxy as well as in more distant (e.g., starburst) 
systems. 

The techniques can be broadly split into four categories. 
The first compares the integrated broad/narrow band pho- 
tometry (e.g., Anders et al. 2004), integrated spectral shape 
(e.g., Smith et al. 2006), or spectral line strengths (e.g., 
Trancho et al. 2007) to a grid of models (age, extinction, 
metallicity) that follow the evolution of the spectral energy 
distribution of a simple stellar population. The second uses 
resolved photometry of individual stars in order to make 
colour-magnitude diagrams (as a proxy for the more physi- 
cal Hertzsprung— Russell diagram) in order to compare with 
theoretical stellar evolutionary isochrones. This technique 
has a long history in clusters within the Milky Way (see, for 



example, Russell 1914) and has recently been extended to 
massive > 10 5 Mq YMCs in nearby galaxies (e.g., Larsen et 
al. 2011). The third category is to use some feature of the 
YMC, such as surface brightness fluctuations, the size of the 
bubble created by the YMC in the surrounding ISM (e.g, 
Whitmore et al. 2011), or the equivalent width of recom- 
bination lines (e.g. Leitherer & Heckman 1995) as a proxy 
for the cluster age. The fourth and final broad category is 
the use of specific spectral features that are due to a unique 
stellar evolutionary phase with known beginning time and 
duration. 

While each technique requires accurate models of stellar 
evolution and spectral synthesis— of which our understand- 
ing and techniques continue to develop— the main difference 
is observational. Techniques requiring only broadband pho- 
tometry are observationally efficient, but spectra provide 
richer sources of information and may warrant the increased 
exposure time based on the goals of the specific project, i.e. 
age, metallicity, morphology, etc. 

An example of the fourth category of age dating tech- 
niques is the use of Wolf-Rayet spectral features (namely 
the 4650A emission "bump") in the integrated spectra of 
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young clusters (e.g., Bastian et al. 2006; Sidoli et al. 2006). 
This type of technique is particularly attractive, as it avoids 
the complexities of modelling a full stellar population, and 
focusses on a single stellar evolutionary phase. Here, we in- 
troduce a new technique along similar lines, based on the 
evolution of massive stars. We apply this technique to a 
photometric catalogue of young massive clusters in M83. 

Red supergiants (RSGs) are young, massive, luminous 
stars evolving from blue supergiants which rapidly exhaust 
their core hydrogen. The tremendous luminosities of RSGs 
peak at ~ 1/j.m with absolute J-band magnitudes of Mj = 
—8 to —11, rivalling the integrated light of globular clusters 
and dwarf galaxies. 

The near-IR flux of RSGs dominates the integrated 
properties of clusters which contain them. For example, a 
single RSG of 25 M at 4000 K emits 25% of the J band 
flux of a 10 s Mq cluster when initially formed. When the 
massive stars in such a cluster begin to evolve to the RSG 
evolutionary phase they emit in excess of 95% of the J band 
light of the full population. In this letter we discuss a new 
technique using near-IR colours as an absolute age indicator 
for young star clusters. This technique utilizes the initial ap- 
pearance of RSGs and their strong effect on the integrated 
properties of their host populations. 

In §2 we describe our observational data, in §3 we 
present the results of our age diagnostic simulations, and 
in §4 we offer brief conclusions and remark on future appli- 
cations of this and similar near-IR techniques. Cluster simu- 
lations in this letter utilize a Saltpeter initial mass function 
with mass boundaries of 0.8-100 M Q 



2 OBSERVATIONS 
2.1 Cluster catalogue 

The catalogue of clusters in M83 was taken from Bastian et 
al. (2011, 2012). We refer the reader to those papers for de- 
tails on the cluster selection, photometry and derivation of 
cluster properties (ages, masses, extinctions, sizes). Here we 
briefly outline the methodology adopted. Cluster candidates 
were selected through visual inspection of Hubble Space Tele- 
scope (HST) Wide Field Camera 3 (WFC3) V-band images, 
from a larger list of candidates selected through the auto- 
mated SExtractor (Bertin & Arnouts 1996) routine, to be 
resolved, centrally concentrated and symmetric in the inner 
regions of the cluster profile. The cluster properties were 
derived through comparison of each clusters U (F336W), 
B (F438W), V (F555W or F547M), Ha (F657N), and I- 
band (F814W) magnitudes with simple stellar population 
models (SSP) using the method described in Adamo et al. 
(2010a). For the SSP models we adopted those of Zachrisson 
et al. (2011) that include nebular emission, and we adopted 
a metallicity of 2.5 times solar (Bresolin & Kennicutt 2002, 
Bresolin et al. 2005). 

We emphasise that only optical and UV colours were 
used in deriving the age of each cluster. Additionally, we only 
included cluster candidates that had inferred photometric 
masses above 5 x 1O 3 M0. This cut was made to limit the ef- 
fects of a stochastically sampled stellar initial mass function 
within the clusters. Catalogue clusters near the minimum 
mass may be affected slightly. In clusters of lower mass, the 



population of massive stars poorly samples the full prob- 
ability density distribution of stellar masses such that the 
appearance of RSGs is skewed to older ages by gaps in the 
cluster population. 

2.2 HST WFC3 near-IR photometry 

In the present work we also make use of the WFC3 Early 
Release Science images of M83 taken in the F110W (J) and 
F160W (H) bands. We carried out aperture photometry on 
these near-IR images, using the cluster catalogue described 
above, with an aperture radius, background radius and an- 
nulus size of 2.5, 4, and 1 pixel, respectively, the same (in 
arcsec) as was used for measuring the optical colours. The ze- 
ropoints were taken from the WFC3 HST Zeropoint website. 
Aperture corrections were derived from bright and isolated 
sources (clusters and individual— likely foreground— stars) . 
We applied aperture corrections of 0.53 and 0.59 magnitudes 
in the F110W and F160W bands. We note that the individ- 
ual aperture correction for each sources varied, depending 
on the extent of which the source was resolved. However, 
the difference between the F110W and F160W bands was 
always 0.06 (with scatter less than 0.03 mag). Finally, we 
only included clusters in the sample discussed here that have 
photometric errors less than 0.1 mag in both near-IR bands. 



3 RESULTS 

3.1 Age vs. near-IR colour 

In Figure [T] we show the F110W-F160W colour vs. the de- 
rived age (from optical colours only) for clusters in the in- 
ner field (top panel) and outer field (bottom panel). In both 
fields we see a bi-modal colour distribution, with clusters 
with ages < 5 Myr having F110W-F160W colours < 0.4, 
and older clusters having redder colours. In particular, we 
note the appearance of a step-like change in colours at 
~ 5 Myr. As will be discussed below, this jump is predicted 
from stellar evolutionary models due to massive stars enter- 
ing the RSG phase, becoming extremely luminous and very 
red. 

The data plotted in Figure [T] include the effects of red- 
dening, which are relatively small in the infrared. Using the 
Cardelli et al. (1989) reddening law and assuming Ry=3.1 
we estimate <5(J-H) ~0.09*Ay. This means that a selective 
increase in Ay by 9 mag is needed for the older clusters 
exactly at the age of the J-H blue vs. red transition to pro- 
duce an 0.8 mag jump in the colors. In principle a YMC 
that is heavily embedded may be extincted enough to appear 
"old" in the J-H index. However, studies of YMCs in merg- 
ing galaxies have shown that the deeply embedded stage is 
relatively short lived, and that nearly all YMCs detected 
in the radio continuum are also detected in the optical, so 
this is unlikely to present a strong bias (Whitmore & Zhang 
2002). The inclusion of emission line (e.g. Br7) photome- 
try in such an analysis would allow such deeply embedded 
clusters, if present, to be detected. On the other hand, the 
comparison with the SSP models and simulated J-H colors 
(see below) is, of course, affected by the reddening. We note 
that the average Ay is about 0.37 and 0.30 for the cluster 
catalogues of the inner and outer fields, respectively. Thus 
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Figure 1. The derived ages of the cluster sample (using optical 
colours) vs. the J — H (F110W— F160W) colour for the inner field 
(top panel) and outer field (bottom panel) fields of M83 (Bastian 
et al. 2011, 2012). These colours have not been dereddened. The 
solid red curve displays our theoretically calculated colour evo- 
lution using Geneva evolutionary tracks and colours calculated 
from Kurucz models for the main sequence and MARCS models 
for RSGs. The other curves are obtained with the GALEV SSP 
code for twice solar (dashed black), solar (red dashed-dotted), 
and half solar (blue dotted) metallicity. Note that the GALEV 
simulations begin at cluster age of 6.6 [log years] while our simu- 
lations end at 7.7 [log years]. The solid and dotted vertical lines 
show the appearance of RSGs and uncertainty as calculated in 
this work. 



Figure[T]is not strongly affected by reddening. The J-H tech- 
nique may be inapplicable to cases of high or strongly varied 
extinction, such as Arp 220 or M82. 



3.2 The effect of RSGs on near-IR cluster colours 

As a massive star evolves into a RSG following the blue su- 
pergiant phase it expands and cools while the bulk of emit- 
ted flux shifts from optical to near-IR wavelengths. These 
objects evolve through the main sequence quickly— within 
10 Myr, for masses above 15 Mq— and through the blue 
supergiant phase even faster, in less than 1 Myr for simi- 
larly massive stars (Meynet & Maeder 2000, Kudritzki et 
al. 2008). 

In Figure [2] we show the dominance of RSGs on the J- 
band flux of massive star clusters as a function of age. These 
simulations are done using Geneva evolutionary tracks with 



Figure 2. The J-band flux contribution of RSGs in a simulated 
super star cluster of 10 5 Mq using Geneva evolution tracks at 
solar metallicity including stellar rotation. Error bars show esti- 
mated scatter from ten simulated clusters using a Salpctcr initial 
mass function. Horizontal lines represent the flux contribution 
from a single 25 Mq 4000 K RSG injected into simulations of 
zero age clusters with masses of 10 4 Mq (blue, dash-dotted), 10 5 
Mq (green, dashed), and 10 6 M (red, dotted). 



rotation at solar metallicity. The work clearly shows that 
the appearance of RSGs marks a significant change in the 
observational properties of the integrated light of the cluster 
in the near-IR. 

The appearance of the first RSG in a young massive 
cluster initializes a hysteresis effect in which near-IR colour 
reddens significantly and remains red as increasing numbers 
of less luminous stars evolve to red evolutionary phases. As 
a single RSG affects the near-IR colours of a massive star 
cluster in this way, the measurement of that change in colour 
represents an absolute age indicator. 

Observations show that 50-60% of massive stars form 
and evolve as binaries (Sana et al. 2012, Sana & Evans 2010). 
The technique presented in this letter requires that the age 
of the first RSGs remain stable (Figure|3J) and that the RSGs 
dominate the near-IR flux of the YMC (Figure [2]). Eldridge 
ct al (2008) argue that, for a reasonable distribution of bi- 
nary separations and mass ratios, the mean lifetime of RSGs 
may be reduced by a factor of 2 to 3 compared with the pre- 
dictions of single-star models. Evolution towards the RSG 
phase is unaffected. The most extreme effect of binarity on 
the stellar population is therefore to remove ~60% of the 
RSGs, to be replaced by blue supergiants. Though this will 
cause the YMC to be fainter in the near-IR, the RSGs' con- 
tribution to the cluster's total near-IR flux decreases by only 
a few percent, and the sharp jump in J-H colour for ages of 
>5Myr remains. Hence, the impact of binarity on our re- 
sults is merely to increase stochasticity effects in low-mass 
clusters. 



3.3 Age dating clusters based on the 
presence/absence of RSGs 

Our simulations localize the colour break in clusters with 
and without RSGS to 5.7 ± 0.8 by weighting multiple tech- 
niques present in the literature to measure the timescale of 
formation for RSGs. 
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Figure 3. Effects of rotation on the mass and age of the first 
RSGs in a simple stellar population. Points represent the first 
appearance of RSGs as a function of initial rotational velocity, 
Vrot (Brott et al. 2011). The solid and dotted black lines show 
the average and scatter in these initial ages excluding the point 
at Vrot = 550 [km/s]. The red dashed line marks the appearance 
of RSGs in Geneva evolution tracks without rotation, and the red 
dash-dotted line Geneva with rotation included at V ro t = 0.4 V c 
(~ 250-350 km/s for the initial population of RSGs). Utilizing 
the three sets of stellar evolution models, the initial appearance 
of RSGs occurs at 5.7 ± 0.8 Myr. 
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Figure 4. Colour evolution of RSGs in temperature and metal- 
licity used in this work. Colours are measured using a new grid of 
MARCS models and TURBOSPECTRUM (Alvarez & Plez 1998, Plez 
2012). 



The uncertainty inherent in the RSG age dating tech- 
nique arises from any process which affects stellar evolution. 
Here we describe the role of stellar rotation, which affects 
the maximum stellar mass which evolves to the RSG phase. 
We utilize the Brott evolutionary models for this discussion 
as they cover a range of metallicities and rotation (Brott et 
al. 2011). Figure [3] shows the evolution of the appearance of 
the first RSG in the Brott evolution tracks as a function of 
initial stellar rotation, as well as the appearance of RSGs in 
two sets of Geneva evolution tracks with and without rota- 
tion at solar metallicity (Meynet & Maeder 2000). 

The Brott models show a slow increase in the time to 
the appearance of a RSG which evolves with a decreasing 
mass of the initial RSGs. This effect is due to an increasing 



size in the hydrogen burning core of fast rotating stars. At 
high rotational velocities, stars evolve as if they were more 
massive, and as such explode in core collapse supernova be- 
fore reaching the RSG phase. At rotation velocities above 
500 km/s this effect becomes significant, lowering the mass 
of the first RSG to ~15M from the consistently 30-40 M Q 
initial RSG population for less extreme rotation velocities. 
Assuming a spread of velocities in a simple stellar popula- 
tion, we drop the 550 km/s case as all stars would not be 
expected to rotate with such extreme velocities and calculate 
an age for the appearance of RSGs to 5 ± 0.7 Myr from the 
Brott models. The Geneva tracks yield slightly longer delays 
before RSG appearance of 5.6 Myr for tracks without rota- 
tion and 6.9 Myr for tracks with rotation. The Geneva group 
models with rotation adopt initial rotation speeds of V ro t~ 
0.4V C . Weighting these three measurements we measure a 
time of RSG appearance of 5.7 ± 0.8 Myr. 

As discussed in §3.2, the effects of these initial popula- 
tions of RSGs on the near-IR cluster colours are expected to 
appear as a significant hysteresis jump allowing for a clear 
differentiation between young "blue" clusters and "red" clus- 
ters in which massive stars have evolved into RSGs and dom- 
inate the near-IR properties. Our cluster catalogue shows 
exactly such evolution in the HST J-H (F110W-F160W) 
colour as plotted in Figure [T] The ages of these clusters 
have been calculated by an independent method utilizing 
optical colours (see §2), and the work in this letter shows 
that a clear age differentiation may be measured using only 
near-IR colours. We model this effect by simulating clusters 
using the Geneva stellar evolution models including rotation. 
Cluster colours are calculated at every time step with syn- 
thetic photometry of cluster stars, utilizing the Kurucz grid 
of spectral energy distributions for main sequence stars and 
a new grid of MARCS models calculated for the critical RSG 
population (Gustafsson et al. 2008). In Figure [4] we plot the 
new grid of colours applied to RSGs in our synthetic cluster 
calculations. This technique successfully recreates the form 
and scale of the observed effect in M83 (Figure [1]). We re- 
cover nearly identical results using the GALEV code to model 
our simple stellar populations over a range of metallicities 
from half to twice solar, and note that the effect is recovered 
in similar passbands using starburst99 (Kotulla et al. 2009, 
Leitherer et al. 2010). 

We note, however, that SSP calculations at metallic- 
ties lower than half solar with different sets of evolutionary 
tracks give conflicting results with regard to the steepness of 
the blue to red (J-H) transition. This is obviously the result 
of a different treatment of mass-loss and rotational mixing. 
A careful observational investigation of stellar clusters in low 
metallicity environments will be worthwhile in this regard. 

Studies such as Sharma et al. (2007) have found a sig- 
nificant spread in the age of the stellar population in young 
associations. Naylor (2009) shows that a systematic differ- 
ence in age determination between young main sequence and 
pre main sequence stars can explain the apparent spread 
in ages. Regardless, our catalogue has been specifically de- 
signed to exclude sources such as the one in Sharma et al. 
(2007). For those true YMCs which remain, no evidence for 
such age spreads in the massive star population have been 
found (Kudryavtseva et al. 2012). 

While near-IR colours represent an efficient method to 
separate young and evolved clusters, spectroscopy in the 
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near-IR offers similar promise to constrain age. In addition, 
spectroscopic observations may be used to extract abun- 
dances of Fe and a-elements of clusters observed at even 
modest resolutions of R~3000 (Davies et al. 2010). This 
is due to the presence of strong and well separated RSG 
spectral features in the near-IR and the dominance of those 
objects in the near-IR (Figure [2]). In two upcoming papers 
we discuss observed and theoretically calculated spectra to 
show the effects of the appearance of RSGs and the possi- 
bility of using such integrated population spectroscopy to 
derive metallicity information (Westmoquette et al. in prep, 
Gazak et al. in prep). 



4 CONCLUSIONS 

In this letter we present a new absolute age indicator for 
young massive star clusters using J-H colour which can be 
used to determine if a massive star cluster is older or younger 
than 5.7 ± 0.8 Myr. Because this technique requires only 
near-IR photometry it can be used well beyond the HST era 
when both space— and ground— based telescopes move to 
the near— IR. The technique can be used with spectroscopic 
data as well and we show this in two upcoming papers. While 
this age dating method is shown to work well in the range 
between half to twice solar metallicities, a careful observa- 
tional investigation of its applicability at significantly lower 
metallicity is suggested, since present evolutionary tracks 
give conflicting results in this low metallicity domain. 
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